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DOTATION 

x,y    Rectangular  coordinates,  x  in  the  longitudinal  direction,  y  in  the 

direction  of  deflection. 
xn ,  y  ,  Rectangular  coordinates,  for  the  vertical  member  of  the  frame. 
Xp,  y?,  Rectangular  coordinates,  for  the  horizontal  member  of  the  frame. 
S  ,  S_,  Shear  force  in  section  1  and  2. 
IL,  M.,  Bending  moment  in  section  1  and  2. 
w      Lateral  distributed  load. 
A      Cross  section  area, 
p      Mass  density. 
E      Young's  modulus. 
I      Second  moment  of  area  of  the  cross  section  about  the  neutral  axis 

through  its  centroid. 
I  ,  I?,  I  for  the  vertical  and  horizontal  member  of  the  frame  respectively, 
to      Natural  frequency  of  vibration  of  the  frame. 
X      Deflection  of  the  vertical  member  at  x  _  L . 


k      A  parameter  equal  to   /- £ 


k.  Corresponding  to  k  for  the  vertical  member,  k  =  ^/iLi£Ai  1 
1  1   y   KL]^ 

k  Corresponding  to  k  for  the  horizontal  member  k  ~  ^  ol^AL  2 

0  Equal  to  k  L  . 

0  Equal  to  k  L  . 
*2  2  2. 


INTRODUCTION 

The  problem  of. computing  natural  frequency  of  vibration  of  structures 
is  of  great  practical  importance  in  the  analysis  and  design  of  structures. 
If  the  frequency  of  the  disturbing  force  is  sufficiently  close  to  one  of  the 
natural  frequencies  of  the  structure,  the  structure  will  undergo  vibrations 
which,  in  the  absence  of  appreciable  damping  effects,  may  become  prohibi- 
tively large. 

The  determination  of  the  natural  frequencies  of  the  lateral  vibrations 
of  beams  has  been  1ixe  subject  of  numerous  papers  and  books.  The  method 
usually  used  in  engineering  applications  is  known  as  the  energy  method,  or 
the  methods  of  Rayleigh  and  Ritz.  A  method  known  as  Transfer  Matrices  has 
grown  out  of  the  well-known  method  of  Holzer  and  Kyklestad.  As  applied  to 
the  determination  of  natural  modes  and  frequencies,  it  is  characterized  by 
the  setting  up  of  a  frequency  determinant  that  follows  as  a  consequence  of 
satisfying  boundary  conditions.  The  computations,  which  are  conveniently 
carried  out  in  matrix  operations,  require  the  insert:./.-  of  trial  values  of 
the  frequency  into  the  transfer  matrices.  Another  powerful  method  of 
analysis  has  been  applied  to  frames  by  Bishop.  Described  as  a  method  of  re- 
ceptanses,  it  deals  with  the  forced  vibrations  of  systems  in  terms  of  the 
receptances  of  its  component  parts.  For  determining  natural  modes  and  fre- 
quencies of  the  structure,  those  frequencies  are  sought  for  which  the  over- 
all receptance  vanishes.  Veletsos  and  Newmark  have  devised  a  method  de- 
scribed as  an  extension  of  Holzer's  method  for  continuous  beams;  it  is  also 
applicable  to  frames  without  sideway.  For  natural  frequency  determination, 


this  method  seeks  those  frequencies  for  which  the  exciting  couple  applied  to 
the  beam  vanishes.  This,  as  in  the  methods  already  mentioned,  requires  re- 
peated trial  computations. 

In  contrast  with  the  above  methods,  the  procedure  described  in  this 
paper  leads  to  the  exact  solution  of  the  Bernoulli-Euler  equation  for  the 
vibration  of  a  rectangular  fixed  end  frame.  The  frequency  equation  is  ob- 
tained from  direct  expansion  of  the  frequency  determinant  that  follows  as  a 
consequence  of  satisfying  boundary  conditions.  In  solving  for  the  eigen- 
frequency  from  the  frequency  equations,  the  1620  digital  computer  is  used. 
The  computer  program  and  results  are  shown  in  the  Appendix  of  this  report. 
It  is  known  that  for  a  structure  of  infinite  degree' of  freedom,  there  is  an 
infinite,  number  of  vibration  modes .  In  this  report,  only  the  natural  fre- 
quencies and  mode  shapes  of  the  first  ten  modes  of  vibration  of  a  simple 
rectangular  fixed  end  frame  are  presented.  However,  with  the  same  proce- 
dure, one  can  find  as  many  modes  as  one  desires.  Rayleigh's  energy  method 
is  also  used  to  find  the  frequency  of  the  first  vibration  mode.  The  result 
is  compared  with  that  of  the  exact  solution,  and  the  error  is  found  to  be 


DERIVATION  OF  BASIC  EQUATION  OF  MOTION 

The  method  of  analysis  which  is  developed  here  is  known  as  the 
Bernoulli -Euler  theory  and  is  based  tipon  the  assumption  that  plane  cross- 
sections  of  a  beam  remain  plane  during  flexure  and  that  the  radius  of  curva- 
ture of  a  bent  beam  is  large  compared  with  the  beam's  depth.  Axial  change  in 
length,  shear  deformation  and  the  rotary  inertia  effect  are  neglected;  only 
the  transverse  bending  deflections  will  be  considered. 

Fig.  1  shows  a  short  element  of  a  beam.  It  is  of  length  $x  and  is 
bounded  by  plane  faces  which  are  perpendicular  to  the  axisj  the  faces  are 
identified  by  the  numbers  1  and  2.  The  forces  and  couples  which  act  upon  the 
element  are  also  shown  in  the  figurej  they  are  the  shear  forces  S_,  S^,  the 
bending  moments  1^,  Mg  and  the  applied  lateral  load  w6x.  If  the  deflexion  of 
the  beam  is  small,  as  the  theory  presupposes,  then  the  inclination  of  the 
beam  element  from  the  unstrained  position  is  also  small.  Under  these  condi- 
tions, the  equation  of  motion  perpendicular  to  the  axis  Ox  of  the  undeflected 
beam  is 


A 


S1  -  S2  +  w&c  =  (pA&c)  ^X  (1) 


where  A  is  the  area  of  cross-section,  p  is  the  mass  density  and  y  is  the  de- 
flexion. If  the  equation  is  divided  by  6x  and  if  the  length  of  the  elemsnt 
is  then  made  indefinitely  short,  the  equation  of  motion  is  obtained  in  the 
form 

£♦--*§  (2) 
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Fig.  1 

The  equation  for  the  rotational  notion  of  the  element  about  Section  1 
may  be  written 

526x  +  M  -  M  +  kw*6x2  -  0  (3) 

where  w*  is  the  mean  value  of  <ju  distributed  on  6x  and  k  is  the  distance  from 
the  center  of  the  distributed  load  co  to  section  2.  When  the  length  of  the 
element  is  made  small,  this  equation  gives 


Now  from  the  simple  bending  theory  we  have  the  relation 


00 


M  =  EI  il  (<Q 

where  E  is  Young* s  modulus  of  the  material  and  I  is  the  second  moment  of  area 
of  the  cross-section  about  the  neutral  axis  through  its  centroid.  If  this  is 
combined  with  equation  (J;)  it  gives 


S  „  -  EI  LZ  (6) 

This  expression  for  S  may  be  substituted  in  equation  (2)  which  then  gives 
the  differential  equation  of  motion  in  the  form 

t?**i^-*t  (7) 

'When  the  applied  lateral  load  w  is  equal  to  zero,  equation  (7)  becomes 

S  &m  .&  (8) 

AP  *J*  ?t2 

For  the  normal  mode  of  vibration,  it  is  reasonable  to  set 

y(t,x)  =  T(t)x(x)  (9) 


and 


a*7     w  ,  «2 


fr  =  Z(x)  — ^  T(t)  (10) 

atr  dt^ 

EI    s\r  EI     d^ 

U  r! " T(t)  ^  T x(x>  <u> 

F  °x  dx 

substitute  equation  (10)  and  equation  (11)  into  equation  (8)  and  rearrange 

^  A2 

\  X(x)     -4,  T(t) 

EI  dx  dt  ,   . 

rP   '  i(x)   =  — W  (12) 

The   separability  of  x  and  t  in  equation  (12)  proves  the  assumption  of  equation 
(9)  is  right.  The  left-hand  member  is  a  function  of  x  alone  and  the  right- 


hand  member  is  a  function  of  t  alone,  the  only  way  for  them  to  he  equal  is 
both  equal  to  a  constant,  say  u)2.  Therefore,  equation  (12)  can  be  reduced  to 
the  solution  of  two  equations.  Namely 

d2        9 

-~  T(t)  +  a)2  T(t)  =  0  (13) 

dt^ 

and 

4xW-^iW-o      •  .  Oh) 

Let 

and  substitute  equation  (lf>)  into  equation  (lit)  to  get 

J* 
■A-  X(x)  -  k^  X(x)  =  0  (16) 

dxa 

Equation  (13)  and  (l6)  are  the  differential  equations  for  normal  vibra- 
tion of  beams. 

The  solution  for  equation  (13)  is 

T(t)  =  F  sincut  +  Fjzosu*  (17) 

or 

T(t)  =  F3sin  (cut  +  e)  (18) 

where  F^,   Fg,  F.  and  9  are  constants  to  be  determined  from  initial  conditions. 


The  natural  frequency  is 

f  .  JtL  (19) 

2tt 

The  solution  for  equation  (l6)  is 

Z(x)  =  G.sinkx  +  G2coskx  +  G-sinhkx  +  G,  coshkx  (20) 

where  G  ,  G  ,  G  and  G,  are  constants  to  be  determined  from  "boundary  condi- 
1  2'     3     h 

tions. 

Equations  (17)  or  (l8)  and  (19)  are  the  basic  equations  for  normal  vi- 
bration of  beams. 


DERIVATION  OF  1HE  FREQUENCY  EQUATION  FOR  NORMAL 
VIBRATION  OF  A  RECTANGULAR  RIGID  FRAME 

A  rectangular  rigid  frame  is  simply  three  straight  beams  joined  at 
right  angles,  therefore,  the  vibration  characteristics  of  a  beam,  as  has 
been  shown  in  the  preceding  paragraphs,  can  also  be  applied  in  frames. 


EI.     l* 
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Fig.  2   Symmetrical  Mode 


Fig.  3   Anti-symmetrical  Mode 


Fig.  2  and  Fig.  3  show  the  two  different  kinds  of  normal  mode  vibrations 
of  a  rectangular  rigid  frame,  namely,  the  synsnetrical  mode  and  the  anti- 
symmetrical  mode.  This  report  will  discuss  these  two  systems  separately. 
However,  in  both  cases  the  method  of  superposition  is  freely  used.  Owing  to 
the  symmetry  of  the  frame  itself,  only  one  of  the  two  vertical  members  will 
be  taken  into  consideration. 

1»   Symmetrical  modes  of  vibration  of  a  rectangular  rigid  frame 

For  normal  modes  of  vibration  of  beams,  the  deflection  curve  is 


y(t,x)   =  T(t)  X(x) 


where 


T(t)  =  F-jSin  (cob,  +  e) 
the  deflexion  curve  for  vertical  bar  a-b,  can  be  written  as 
X.  (x)  »  A,  sin3c,x,  +  A^oskjX,  +  A^sinhkjX,  +  A>  coshkjX, 


(9) 


(18) 


(1-1) 


H^-^—^ 


H, 


Vertical  bar  vertical  bar 

Fig.  h 

the  deflexion  curve  for  the  horizontal  bar  b-c  can  be  written  as 

Z2(x2)  =  Ej^ink^  +  B2cosk2x2  +  BosinhkgXg  +  Bi  coshk^       (1-2) 

where  x_,  x  ,  y  ,  y  ,  are  coordinates  as  shown  in  Fig.  k,   k. ,  k?,  are  the 
corresponding  to  k  for  bar  a-b,  b-c,  in  equation  (20).  A_,  Ap,  A.,  A,  and 
B  ,  B2,  B  ,  B,  are  constants  to  be  determined  by  boundary  conditions. 
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Boundary  Conditions 

(1)  For  the  vertical  bar  a-b. 

(a)  X^) 

=  0 
Xl=° 

(b)  IjCx^ 

=  0 

x1  =  0 

(c)  X^) 

=  0 

A2  -  -  \ 


*l"-^J 


(1-3) 

U-M 


Ao  (sinhk.  L.  -  sink.  L.  )  +  A.  (coshk.  L  -  cosk.  L, )  =  0 


coshk,  L,  -  cosk-,I», 


(1-5) 


3    sinbk-jL-,  -  sink-,L,  *~h 
substitute  equations  (1-3),  (l-U)  and  (1-5)  into  equation  (l-l)  and  simplify 


i-  coshk-iL,  -  cosk-iL-i  -i 

X(x  )  =  A  coshk  x  -cosk  x  +  i~± ±-i  (sink  x  -sinhk  x  ) 

1  1    UL    11    11   sinhk^  -  sink^     11     1  1  -' 


(1-6) 


differentiate  equation  (1-6)  and  simplify 


r  coshk,  I»  -cosk-,  L,  -i 

X'(x  )  =  k  A.  sinhk  x  +  sink  x  +  ±~± £-£.  (cosk  x  -sinhk  x  ; 

11  1  UL  11  11       sinhk  L  -sink  L  11  1  1  - 

(1-7) 

differentiate  equation  (1-7)  and  simplify 

2     r  coshk,  L-,   -  cosk^L,  -, 

(1-8) 

(2)    For  the  horizontal  bar  b-c 


(a)     X,^) 


=  0 


x2  =  0 


B2  =   -  Bk 


(1-9) 
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(b)     X^) 


=   0 


x2  =  L2 


substitute  equation  (1-9)  into  equation  (1-2) 

B1sink2L2  -  B^ooskgLg  +  B^sinhkpl^  +  B^coshkgl^  =  0 


(1-10) 


(c)  X2(x2) 


=  0 


=  X^(x2) 


Xg  -L, 


Combine  boundary  condition  (c)  and  equation  (l-9) 

-B-sink  1^  +  B, coskgL  +  B^sinhk  Lg.  +  B,  coshlc^Lg  =  2B,     '  (1-11) 
solve  equation  (1-10)  and  equation  (l-ll)  for  B1  and  B,  in  terms  of  B,  , 

(1-12) 


B    cosk2L2  ^r 
1=  " 


sink2L2    ^ 


1  -  coshkJLp 
B3  =   sinhk0L0   BU 


(1-13) 


^2 


substitute  equation  (1-9),  (1-12)  and  (1-13)  into  equation  (1-2) 

.     .  ,coskpLo-l  1-coshkpLo 

X2(x2>  "  V  sm.CL    sink^-cosk^  +    w  j       (sinh^x^oshk^) 


'2^2 


Differentiate  equation  (1-lU)  and  simplify 


(l-H;) 


coskpLp-l  l-coshkplo 

X2(x2)  -  k2Blt(^SE-I_  cosk^+sink^  +-EIEH^-    coshk^sinhkgX,,) 

(1-15) 

differentiate  equation  (l-13>)  and  simplify 


12 


2     f     cosk2L2-l 


l-coshk  L 


X^(x2)  =  k^-     sil±2L2     sink^cosk^  +     six^h     sinhk^coshk^) 

(1-16) 


Compatibility  conditions. 


(1)  ■  X1(x1) 


X2(x2) 


x2  =0 


x1  =  L1 
equating  equation  (1-7)  and  equation  (l-l5)  and  simplify 


k„A 


2(coslc_L1coshk1L1  -  1)  cosk2l<2  -  1      1  -  coshkgLg 


l\      sinhk^  -  sink^      "    2  k  K     sink^      +       sinhk^ 


(2)  EE^XjCx^ 


x  =  L 
11 


-  k2Bl- 
=  EI  X"(xJ 


)-  0 


(1-17) 


"2  2*  2' 


*2.0 


substituted  from  equation  (1-8)  and  equation  (1-16)  and  simplify 
2(cosk,Lsiiihk.1L,  -  sink-,L,coshk,L-, ) 


sinhk  L.  -  sink1L1 


Wi  -  W2 


Irt  =  0 


(1-18) 


solve  A,,  B.  from  equations  (l-17)  and  (1-18).  Both  equations  (1-17)  and 
(1-18)  are  homogeneous  equations.  In  order  to  have  a  non- trivial  solution, 
the  determinant  of  these  equations  must  equal  to  sero. 


A- 


2 (cosk^L-jCoshk-jIn  -  l) 
sinhk  L    -  sink  L  1 


cosk2li2  -  1      1  -  coshkpL/j 
~2(      SISkgLg    "  ""  sinhkj^     ' 


2(cosh^L^sinhk.jL1  -  sink-^L^coshknLi )     2 


sinhk  L  -  sink  L 


-  Hh 


simplify  and  rearrange 
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sir3k2L2      +       sinhkgLg       =    2 


coskgl^  -  1      1  -  coshk2L2  k2I2  cosk^LjCoshk-jI^  -  1 

k_I_     cosk-Lsinhk-L,   -  sin^L-cosbk-L, 

(1-19) 


Let  k.1-  =  0    and  k2L2  =  fl2,   substitute  into  equation  (1-19) 

cos0?  -  1   1  -  cosh02    k  I     cos0  cosh0  -  1 

siii0   +   sinhp?   =   k-i_  cosP-sinhp.  -  sinjZLcosh^. 

(1-20) 

equation  (1-19)  or  equation  (1-20)  is  the  frequency  equation  for  symmetrical 

modes  of  vibration  of  rectangular  rigid  frame. 

2.   Anti-symmetrical  modes  of  vibration  of  a  rectangular  rigid  frame 
For  normal  modes  of  vibration  of  beams,  the  deflection  curve  is 

y(t,x)  «  T(t)  X(x)  (9) 

■where 

T(t)  =  F  sin  (a*  +  0)  (18) 


Hi   U 


ace. 


I, 


-I 


f 


///// 


V. 


Fig.  5 


Ik 


The  deflection  shape  for  the. vertical  bar  a-b  can  be  written  as 

X^(x  )  =  C  sink  x  4.  c  cosk  x  +  C  sinhk^  +  C.coshk^   (2-1) 

The  deflection  shape  for  the  horizontal  bar  b-c  can  be  written  as 

XgCXg)  =  D^inkgXg  +  DgCosk^  +  D-sinhk^  +  D^coshk^   (2-2) 

where  X,,  X2,  7-,,  y2,  are  coordinates  as  shown  in  Fig.  J>,  k^,  k2  are  the 
corresponding  to  k  for  bar  a-b,  b-c  in  equation  (20).  C^,  Cg,  C\_,  Ci  ,  and 
D_,  D«,  D^,  Di  ,  are  constants  to  be  determined  by  boundary  conditions. 


Boundary  conditions 

(l)  For  the  vertical  bar  a-b 

! 

(a)  X^) 

=  0 

(b)  X^) 

=  0 

C2  =  -  C^         (2-3) 


C1  =  -  C3         (2-li) 


substitute  equation  (2-3)  and  equation  (2-li)  into  equation  (2-1)  and  simplify 
£,00  =  C^sinhk-si  -  sink^)  +  C,  (coshk^  -  cosk^)       (2-5) 

differentiate  and  simplify 

X«(x  )  =  IslC  (coshk-JC-  -  cosk  r  )  +  k  C,  (sinhts:    +  sink^O         (2-6) 

differentiate  equation  (2-6)  and  simplify 

X"CO  =  kPc  (sinhk^  +  sink^)  +  k^C,  (coshk^  +  cosk^)         (2-7) 

differentiate  equation  (2-7)  and  simplify 

X^K^)  =  k^C  (coshk^H-  cosk^)  +  k^C,  (sinhlc^  -  sink^)    (2-8) 
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(2)  For  the  horizontal  bar 

=  0 


(a)  X2(x2) 


(b)  X2(x2) 


x2  =  0 


V-Du 


(2-9) 


Xg  -  0 

combine  equation  (2-9)  and  equation  (2-2) 

-Dl sinkgLg  -  D,  coskgl^  +  D^sinhkgLg  +  D^coshkgLg  =  0 

=  -X»(xJ 


(2-10) 


(c)    X»(x2) 


x2  =  0 


2V  2' 


*2    =L2 


-ILsinkgLg  +  D,  coskgLg  +  D-sinhkgLg  +  D.  coshkgLg  =  -2D^         (2-11) 
solve  equation  (2-10)  and  equation  (2-11)  for  Dp  D~,  in  terms  of  D^ 


1  +  coshkpL,, 
D3  =  sinhk2L~  \ 

1+  cosk2L2 


(2-12) 


(2-13) 


substituting  equations  (2-9),   (2-12),   (2-13)  ixito  equation  (2-2)  and  simplify 


.1  +  cosk-Lg 


X2(x2)    =Bh[(  )  sink^  -  cosk^ 


-  (" 


2  2 

1  +  coshk  L2 


sinhkgl^     5  siahkgXg  +   coshk^] 


(2-lit) 


differentiate  equation  (2-liO  and  simplify 

1  +  cosh2L2 
X«(x2)  =  k2Bh[(     s±^h      )  cosk2x2  +  sinlc^ 


1  +  cos 


tokgLg 


-  (    sinhk  I )  coshkgJ^  +  sinhk^  (2-15) 


16 

differentiate  equation  (2-15)  and  simplify 

1  +  cosk2I2 

Xgfcfe)  -  kgD^-  (     siniC2L2     )  sinkgXg  +  cosk^ 

1  +  coshkgl^  _ 

"  ( — sinhkL ^  sinhk^  +  coshk2X2  '  (2-l6) 

differentiate  equation  (2-16)  and  simplify 

1  +  cosk  L 
T£x2)  =k2Du[-  (     sil&2^  2)  coek^  -  .sixfc^ 


1  +  coshk  L-  -, 

"  <     rtdft^     >  °°Shk2X2  +  sMV2  J  (2"17) 


Compatibility  equations 

(1)  To  satisfy  the  continuity  condition  at  joint  b 

.     \*i  -h  V  ° 

substituted  from  equations  (2-6)  and  (2-5)  and  simplify 

C  k-Xcoshk-L  -  coslcL  )  +  C,k-,(sinhk  L  +  sinlcL.) 

1  +  coskpLp   1  +  coshkpl^ 
"  k2Dit  (  sinkgLg      sinhk^  *  =  ° 


(2)  To  satisfy  the  static  equilibrium  condition  at  joint  b 


substituted  from  equation  (2-7)  and  (2-l6)  and  simplify 


(2-18) 


17 


C  I  k^sinhk  L  +  sinl^)  +  (coshlt^  +  cosk^)  C^kJ  -  2\J2^  =  ° 

(2-19) 

(3)  To  satisfy  the  dynamic  equilibi-ium  condition 

The  force  to  produce  the  horizontal  motion  of  the  horizontal  bar  is  the 

shear  force  at  the  top  ends  of  the  vertical  bars. 

total  shear  force  =  mass  of  horizontal  bar  times  ace.  of 
horizontal  bar. 

3  2 

2EL.  l~-~  y(t,  x. )       1  -  -  (Ap)  L  S-W  y(tp^) 

from  equation  (9) 

7(t,^.)  =  T(t)  X^)  (2-20) 

3 
-*-  y(t^)  =  T(t)  XJC^)  (2-21) 

Sx 

^2  *2 

^Ty(t,xL)  =  J^  T(t)  2^)  (2-22) 

since 

T(t)  =  F-,sin  cut  +  P2cos  cut  (17) 

£-*  T(t)  -  -  U)2T(t)  (2-23) 

dt 

substitute  equations  (2-20),  (2-21),  2-22)  and  (2-23)  into  the  dynamic 

equilibrium  equation  and  simplify 
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C     [ELLl^(coshk1LL+    cOsk^)  +  1/2  (Ap^I^teinhk^  -  sinlCjL^J 

C,    [m^CsinhkjI^  -  sinkjl^)  +  1/2  (Ap)2L2u>  (coshk^  -  cosk^)]  „    0 

(2-210 

solve  C  ,  C,  ,  D,  from  equations  (2-l8),  (2-19)  and  (2-2U).  Both  equations 
(2-18),  (2-19),  and  (2-2k)  are  homogeneous  equations.  In  order  to  have  a 
non-trivial' solution,  the  determinant  of  these  equations  must  equal  zero. 

ki  (eoshk-L.  -  coslc-L.)         kl  (ginlikJL,  +  sink^) 

1   cosk0L0   1   coshk0L0 

k  ( LL LL\ 

V     sink2L2  sinhk2L2     ; 


A  = 


I  k?(sirihk  L     +  sink  L _)  3^  (coshfc^  +  cosk^)  21^ 


EI  k^coshlcl^   +  coslcO       +        ELk_  (sirihkjl^  -  sinkjL^  + 

(Ap)2         2  Up) 

— s —  L2cu  (sinhk-L,    -  sinkJLJ       — __£  L  t«  (coshk  L     -  cos3c  L  ) 


expanding  the  determinant  and  simplify 

-,     1  +   cosk^Lj       1  +  coshk9Lp  0 

*  ■  I  (    **ft      "  — iSEjf")  I1k1„2(A0)2L2(2Stak1l1coShk1l1 

-  sinhk_L_cosk  L ) 

2  ,    1  +  cosk  L  ■      1  +  coshk  L 
+   "ft?     s^24       -  -inhkiLp)   ("2  "  Zcoshk^cosk^) 

+   2EI2k2I  k_(2coshk  L1sink1L1  +  2cosk1L1sinhk1L1) 
-    I2k2(Ap)2L2u?(2  -  coshk^cosk^)  =  0 
rearrange  and  simplify,  and  let  0   =   k,L,,  0   «   k^ 
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I_k_  cos02sinh02  -  sin02cosh02.  +  sinh.0  -   sin0 
Igk^  2sinppSinh02 

(Ap)2L2 

rrr%-|-- 0-^(1  -  cosh0 cos01)  -  2(cosh01sin01 +  cos01sinh01) 

"Ta^ <*-*> 

(Ao")  L  '  &l  (s^fti003*1^  ~  sinh0.cos01)  -2(1+  cosh01cos01 ) 

equation  (2-2f>)  is  the  frequency  equation  for  anti-symmetrical  modes  of  vibra- 
tion of  a  rectangular  rigid  frame. 
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NUMERICAL  ILLUSTRATIVE  EXAMPLE 


For  a  rectangular  rigid  frame  as  shown  in" Fig.  6,  find  the  natural 
frequencies  and  the  shapes  of  vibration  of  the  first  eight  normal  modes. 


K 


U 


HA/% 

4BL, 


(API 


mi 


(An 

EI. 


"*" 


L1=L2 


I. 

•L2 


TaTJJ 


=  o.£o 


Fig.  6 


from  the  data  given 
Jlkl   V2 


2  2 


(Ap)2I2 

and  11*555;  ■  2-o 

s  frequency  equation  of  anti-symmetrical  modes,  equation  (2-25),  takes 

the  form 

krz   cos0  sinh0rt  -  sin0  cosh0    +  sinh0    -  sin0 
v  <Z  2  2  2  £  c.  < 

8  2sin0«sinh09 


0  (l  -  cosh0  cos0  )  -  (cosh01sin01  +  0030^5^110.^ 
0  (sin£i  coshp    -  sinh0  cos^)  -  (1   +  cosh^cos^J 


(3-D 


and  the  frequency  equation  of  symmetrical  modes,  equation  (1-20)  takes  the 


:Jc_-_ 
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C03^2  -  1      1  -  cosh^        g  cosj^coshj^  -  1 

•      sin02      +      einh02       "  17^    cosft  sdUihi  -  sin0  cosh0  '3-2) 

Now,  the  problem  becomes  one  of  finding  the  eigenvalue  u>  to  satisfy  equation 
(3-1)  or  equation  (3-2).     It  can  be  solved  by  trial-and-error  method  with  the 
help  of  the  digital  computer. 
In  this  problem,  note  that 

hm  Tzh  (3-3) 


and  let 


II 


0X(1  -  ooshf^cos^)  -  (coshjijSlnJZ^  +  cos^sinh^) 

p^sinp  coshp  -  sinhJJ  oos0  )   -  (1  +  cosh?  oosp  )     &m) 

V|  oos02sinh02  -  sin02cosh02  +  sinh02  -  sin02 
"  B  2sin02sinh02  "  &m^ 

"o   CO30  sinhp  -  sinp  coshJJ  '3-6 ) 

cos02  -1    1  -  cosh02 
12  "  "~iin^ +   sinh02  (3"7> 

and  then,  set  up  Forego  program  to  find  the  value  of  II,  12,  II  and  T2  for 
^  ranging  from  0,0  to  12.0  with  each  0.1  increment,  the  results  are  shown  in 
Appendix  1.  With  this  results,  curves  for  II,  12  against  0L  (Fig.  7)  and  for 
II,  T2  against  0^   (Fig.  8)  can  be  plotted.  So  the  0  values  for  which 
n  -  12,  and  11  -  T2  can  be  located. from  the  figures.  It  should  be  noted  that 
the  0    value  thus  found  will  be  accurate  only  for  the  first  decimal  place. 
For  example,  the  0    for  the  first  anti-symmetric  node  can  be  between  1.7  and 
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-s.< 


-&.» 


-7c- 


S.c 


■isV 


Fig.   8     Symmetr 


Leal   modes 
Y2    =  horizontal   bar 
Yl   =  vertical   bar 


2U 

1.8.  If  two  decimal  place  accuracy  is  desired,  the  same  program  can  be  used 
to  find  the  value  of  II,  X2,  XI,  Y2,  but  for  01   ranging  from  1.71  to  1.80 
with  each  0.01  increment.  Repeating  this  process,  one  can  find  as  much 
accuracy  as  one  wants.  This  report,  however,  has  reached  only  four  decimal 
place  accuracy.  From  the  relation  fL   ■  k_L  and 


u>  (pA)n 


■s: 


kl 


(15) 


solve  for 


*&* 


(3-8) 


Table  1  shows  the  01  and  ^  values  for  the  first  ten  normal  modes  of  vibra- 
tion of  a  rectangular  rigid  frame  as  shown  in  Fig.  6.  The  value  of  07  can 
determine  the  natural  frequency  u>,  as  expressed  in  equation  (3-8). 

Table  1.   7alues  of  01   and  ^  for  the  first  ten  normal  modes 


modes 

1 

2 

3 

U 

5 

6 

7 

8 

9 

10 

h 

Anti- 

symra. 

1.6775 

U.7187 

7.3196 

8.291U 

10.9617 

Symm. 

3.8063 

U.8888 

7.7136 

10.5998 

H.75U6 

A 

2.81U 

1U.U88 

22.255 

23.900 

53.576 

59.500 

68.7U7 

112.356 

120.159 

138.171 



Now  the  next  step  is  to  find  out  the  deflection  shape  of  the  rectangular 
rigid  frame  corresponding  to  each  normal  mode  of  vibration.  This  is  equiva- 
lent to  solving  the  constants  ^  —  A^,  ^  —  B^  in  equations  (l-l),  (1-2), 


* 
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ci  —  CV  Di  —  \  ^  ec?'uations  C2"1)?  (2-2).' 

For  the  constants  A-, A^,  B, B*  In  equations  (l-l)  and  (1-2) 

first  solve  equation  (l-lS)  for  Ai  in  terms  of  Bk 


sinh01  -  sin^ 


a  c~  . i  R  (3-9) 

4      E7f    (cos01sinh01  -  sinj^coshjy     U 


by  equation  (1-5) 


cosh0L  -  cosjl 

S  =  "  sinh01  -  sin^  AU 


by  equation  (1-3) 


by  equation  (l-U) 

\---s 

by  equation  (1-9) 


by  equation  (1-12) 


by  equation  (1-13) 


,cos02  -  1 
%=   sinp2  Bh 


(3-10) 


A0  =,  -  A.  (3-11) 

2    h 


(3-12) 


B  -  -B,  (3-13) 

2     4 


1  -  cosb.0 


(3-15) 


How  let 

\  =  1  (3-16) 

from  equation  (3-9)  to  equation  (3-l6),  the  constants  A-^ Al,  B^  - —  Br 
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of  equations  (l-l)  and  (1-2)  can  be  solved  and  the  characteristic  shapes  of 
vibration  of  the  rigid  frame  for  symmetrical  modes  are  Imam.  These  values 
are  listed  in  Table  2. 

Table  2.  Constants  for  characteristic  shapes  of  vibration  of  rectangular 
rigid  frame  for  symmetrical  modes* 

Mode    01              A]_      A2      A3       A^  B±  B,,  ^  ^ 

2   3.8063  -10.2697  10.18U8  10.2697  -10.18U8  33.82U7  -1  -0.9217  1 

k       U.8888  l.hh60    -1.U712  -1.UU60   1.1*712  1.8990  -1  -0.9677  1 

6   7.7136   -1.9778   1.9763   1.9778  -1.9763  -0.1019  -1  -0.9970  1 

8   10.5998   3.1309  -3.1310  -3.1309  .'3.1310  -3.821-9  -1  -0.9997  1 

10   11.7&6   1.1896  -I.I896  -1.1896   I.1896  U-27^5  -1  -0.9999  1 

For  the  constants  C^^ G»  ,  E^ D,  in  equations  (2-1 )  and  (2-2)  first 

solve  equation  (2-18)  and  (2-19)  for  C_  and  Ci  in  terms  of  Q  • 

1   1  +  cos02   1  +  cosh02  ft 

im   <~ siE^ '   sinh0o  )(cos^i  +  cosh^)  -  •^(sinh01  +  sinj^) 

(cosh  0-   -  cos  0.)  -  (sinh0.,  +  sin01) 

(3-17) 
g                1               1  +  cos02   1  +  cosh02 
— (coshj^  -  cosjy  -  j^r  (sinh^  +  sii^)  (  sin^ "sinh02  '  } 

cij.  =  2       2  2 

(cosh  0,  -  cos  0,  )  -  (sinh0.  +  sin0, ) 

(3-18) 
by  equation  (2-3) 

C2  =  -  C^  (3-19) 
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by  equation  (2-1+) 


C1  =  -  C  (3-20) 


By  equation  (2-9) 


D2  -  -  D^  (3-21) 


by  equation  (2-12) 


1+  cosh^L 
D3  =  "  -SHS^  Di*  (3"22) 


by  equation  (2-13) 

1+  cos02  . 


and  let 


°i =       siigL     ^  <>23> 


1^-1  (3-21;) 


from  equations  (3-17)  to  (3-21;),  the  constants  C  C,  ,  D  C,  of 

equations  (2-1)  and  (2-2)  can  be  solved  and  the  characteristic  shapes  of  vi- 
bration of  the  rectangular  rigid  frame  for  anti-symmetrical  modes  are  known. 
These  values  are  listed  in  Table  3. 

Table  3.  Constants  for  characteristic  shapes  of  vibration  of  rectangular 
rigid  frame  for  anti-symmetrical  modes 

Mode    ^      ^       ^      A3      \   _    ^  B2  B3    ^ 

1   1.6775   4*.ll*5l   3.1*U35   l*.ll*51   -3.1*1*3$   1.171*5  -1  -1.61*55  1 

3   3.3992   3.3992  -3.2*385  -3.3992   3.1*385  -0.1*381*  -1  -1.0386  1 

5   7.3196  -11.3996  11.3971*  11.3996  -11.3971*  -15.591*1  -1  -1.001*3  1 

7    8.2911*   -2.2955   2.291*2   2.2955   -2.291*2   2.7869  -1  -1.0019  1 

9   10.8917   3.2087  -3.2088  -3.2087   3.2088   0.1039  -1  -1.0002  1 
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with  these  constants  known,  the  characteristic  shapes  of  the  vibrating  frame 
can  he  computed  from  equations  (l-l),  (1-2)  for  symmetrical  modes  and  from 
equations  (2-1),  (2-2)  for  anti-symmetrical  modes  as  listed  in  the  following 
tables. 

For  symmetrical  modes,  the  characteristic  shape  equations  are: 
the  vertical  bar 
X_(x.^ )  =  JLsinl^XL  +  A^coslc-X.  +  A-sinhk^x,  +  Aicoshk-x.         (3-25) 

the  horizontal  bar 
XgCxg)^  RjSinkgXg  +  BgCoskgXg  +  B^sinhk^  +  A,coshk2x2         (3-26) 

For  anti-symmetrical  modes,  "the  characteristic  shape  equations  are: 
the  vertical  bar 
X_(x,)  =  C  sink  x^  +  C.cosk-X  +  C  sinhk  x+  C.coshk  x.         (3-27) 

the  horizontal  bar 

XpCxp^  D-sinl^Xp  +  DpCoskgXg  +  D-sinhlo^  +  DtcoshkpXg         (3-28) 

Table  4  to  Table  13  contain  the  corresponding  ]L  (x.  )  and  !„  (x_ )  values  for 
equation  (3-25)  and  equation  (3-26)  in  tabular  form  for  the  first  five 
symmetrical  modes.  While  Fig.  9  to  Fig.  13  are  the  corresponding  character-, 
istic  curves .  Table  lb  to  Table  23  contain  the  corresponding  XL.  (x,  )  and 
X2(x.)  values  for  equation  (3-27)  and  equation  (3-28)  in  tabular  form  for  the 
first  five  anti-symmetric  modes.  While  Fig.  Ik   to  Fig.  18  contain  the  corre- 
sponding characteristic  curves. 
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First  symmetrical  mode  —  0  =  3. 8053,  0g  =  3.2007 

Table  h.     Calculation  for  equation  (3-25),  fly.  =  3.8053 

12  3  U  •   5         6  7 

k^xi           x/L  A-^inkj^  AgCosk^  A  sinhk^  A^cosh!^^  3+l*+5+6 

0.3806  0.100  -3.8152  9.htt6  k.0Ck2          -10.9313  -1.2867 

0.7613  0.200  -7.0861  7.3728  8.5957          -13.2820  -li.3996 

1.11*19  0.300  -9.3393  it.2359  li*.l*l*7l*         -17.5790  -8.2305 

1.5225  0.1*00  -10.2571;  0.1*919  22.1*167          -2l*J»527  -11.8015 

1.9000  0.500  -9.7182  -3.2927  33.5631*          -3U.8086  -ll*.256l 

2.3000  0.605  -7.6581  -6.7861  50.7015         -51.3029  -15. 01*56 

2.7000  0.710  4i.3893  -9.2801  76.0605          -76.1161  -13.6530 

3.0000  0.798  -1.UU91  -10.0830  102.8808  -102.5375  -11.1888 

3.1*000  0.89k  2.6239  -9.81*67  153.6902  -152.7588  -6.29li* 

3.8063  1.000  6.331*1*  -8.0165  230.8700  -229.1380  -0.0001 
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Table  5.  Calculation  for  equation  (3-26),  02  =  3.2007 


1 

2 

3 

k 

S 

6 

7 

k2*2 

x/L 

B]Lsink2x2 

Bgcosk^ 

B-siohkpXp 

B.  coshk2x2 

3+li+5+6 

0.3201 

0.100 

10.61*1*6 

-0.9192 

-0.3001 

1.0517 

10.1*1*70 

0.61*01 

0.200 

20.2035 

-0.8020 

-0.6311 

1.2120 

19.9821* 

0.9602 

0.300 

27.7126 

-0.573U 

-1.0271* 

1.1*975 

27.6093 

1.2803 

0.1*00 

32.1*071* 

-0.2861* 

-1.5299 

1.9378 

32.5288 

1.600k 

o.5oo 

33.8112 

0.0296 

-2.1905 

2.578U 

3U.2287 

1.9201; 

0.600 

331.7781* 

0.3U25 

-3.O772 

3.U851 

32.5287 

2.2000 

0.687 

27.3173 

0.5885 

-It. 1081 

U.5679 

28.395U 

2.6000 

0.813 

17.1*366 

0.8569 

-6.1705 

6.7690 

18.8920 

2.9000 

0.901* 

8.0909 

0.9710 

-8.3502 

9.111*6 

9.8263 

3.2007 

1.000 

-1.9991 

0.9983 

-11.291*7 

12.2952 

0.0X0 
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Second  symmetrical  node  ~  0.  -  1*.8888,     &,  =  k.1110 

Table  6.     Calculation  for  equation  (3-25),     0^  =  li.8888 

1  2  3  U  5  6                     7 

1C..3C.  X./L  A-sirik^  ApCosk-jX,  A-sinhk-x,  AiCoshk-x.  3+l*+5+6 

0.1*889  0.100  0,6792  -1.2989  -0.735U  1.6505  0.295U 

0.9778  0.200  1.1992  -0.8221  -1.6502  2.232k  0.9593 

1.1*667  0.300  1.1*382  -0.1529  -2.9673  3.3585  1.6765 

1.9556  0.1*00  1.31*03  0.5523  -5.0099  5.3031*  2.1861 

2.1*000  0.1*91  0.9768  1.081*9  -7.901*1  8.1753  2.3329 

2.9000  0.593  0.31*59  1.1*285  -13.1002  13.1i09l*  2.0336 

3.1*000  0.695  -0.3659  0.1*221*  -21.61*00  22.0661  1.1*790 

l*.5O00  0.920  -1.1*135  0.3101  -65.071*3  66.221*7  0.01*70 

1*.7000  0.960  -1.1*1*58  0.0182  -79.1*852  8O.8838  -0.0290 

li.8888  1.000  -1.1*236  -0.2611  -96.001*7  97.6889  -0.O0Q5 
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Table  7.  Calculation  for  equation  (3-26),  0  =  U.1110 


1 

2 

3 

k 

5 

6 

7 

2  2 

XgA 

Blsil3k2X2 

B2cosk2x2 

B  sirihk^ 

B,  coshk-X- 
4            2  2 

3+li+5+6 

o.Um 

0.100 

0.7583 

-0.9167 

-0.1;091 

1.0857 

0.5187 

0.8222 

0.200 

1.3912 

-0.6806 

-0.8883 

1.3575 

1.1798 

1.2333 

0.300 

1.7919 

-0.3311 

-1.5199 

1.8619 

1.8028 

l.Qilik 

o.Uoo 

1.8939 

0.0735 

-2.1019 

2.6855 

2.2ltL0 

2.1000 

0.512 

1.6202 

0.50U8 

-3.8920 

li.lUi3 

2.3773 

2.5000 

0.608 

1.1366 

0.8011 

-5.851*8 

6.1323 

2.2152 

2.9X0 

0.706 

O.U5U2 

0.9710 

-8.7670 

9.11U6 

1.7728 

3.3000 

0.803 

-0.2995 

0.9875 

-13.1006 

13.57W 

1.1622 

3.7000 

0.901 

-1.0061 

0.81-81 

-19.5585 

20.2360 

0.5195 

U.mo 

1.000 

-1.5660 

0.5658 

-29.5107 

30.5120 

0.0001 
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Third  symmetrical  mode  —  0  =7.7136,  02   =6.1863 


Table  8.  Calculation  for  equation  (3-25),  &  =  7.7136 


1 

2 

3 

1* 

5 

6 

7 

■ft 

■iA 

A  sink  x: 

A2COsklXl 

A  sinhk  xl 

AjCoshlcXw 

3+it+5<-6 

0.771U 

0.100 

-1.3787 

1.U168 

1.6815 

-2.59bl 

-0.871*5 

1.5U27 

0.200 

-1.9770 

0.0555 

1*.1*139 

-li.8330 

-2.31*06 

2.3000 

0.298 

-1.1*71*8 

-1.3168 

9.761ib 

-9.9S$0 

-2.9822 

3.1000 

0.ii02 

-0.0823 

-1.971*5 

21.9075 

-21.9791* 

-2.1291 

3.9000 

O.U96 

1.3603 

-l.i*3l*6 

•  1*8.831*1 

-1*8.8369 

-0.0771 

b.6000 

0.596 

1.9653 

-0.2217 

98.3700 

-98.3151* 

1.7982 

5.l*ooo 

0.700 

1.5281* 

1.25UU 

218.91*1*2 

-218.7873 

3.2097 

6.2000 

o.8ob 

O.16U2 

1.969b 

1*87.2775 

-1*86.9119 

2.1*992 

6.9000 

0.895 

-1.11*1*0    • 

1.6121 

981.2596 

-980.5173 

1.2101* 

7.7136 

1.000 

-1.9581* 

0.2765 

2213.7397 

-2212.0588 

-0.0010 

3U 


Table  9.  Calculation  for  equation  (3-26),     02  =  6.I4.863 

1  2  3  h 

ICjX.  Xj/L  B^-sinl^Xp  BpCoskpXU 

O.6U86  0.100  -O.C616  -0.7969 

1.2973  0.200  -0.0981  -0.2701 

1.91^9  0.300  -0.09ii8  O.366U 

2.6000  O.liOl  -0.0525  0.3569 

3.2000  0.1*93  0.0060  0.9983 

3.9000  0.601  0.0701  0.7259 

li.5000  0.693  0.0996  0.2103 

5.2000  0.801  0.0900  -0.1*685 

5.8000  0.89k  0.0it73  -0.8855 

6.1863  1.000  -0.02C6  -0.9791; 


5 

6 

7 

B-sinhkpXp 

B»  coshk^ 

3+U+5+6 

-0.6929 

1.2178 

-0.3336 

-1.6880 

1.9663 

-0.0899 

-3.U182 

3.571U 

0.1t302 

•-6.67li6 

6.7690 

O.8988 

-12.2092 

12.2866 

1.0817 

-21;.  6170 

2li.7113 

0.8903 

-lUi.8680 

ii5.oiia 

0.1565 

-90.3615 

90.6389 

-0.1011 

-l6ii.6529 

165.1513 

-0.3398 

-327.C608 

328.032U 

-O.OOOlt 

3p 


Fourth  symmetrical  mode  —  01  =   10.5998,  02  =  8,9133 
Table  10.  Calculation  for  equation  (3-25),  0,  =  10.5998 


1 

2 

3 

h 

5 

6 

7 

to 

Xj/L 

A^inlc^ 

A2cosklxl 

A,sinhk1x1 

AiCoshkLX, 

3+li+*-6 

1.0600 

0.100 

2.731U 

-1.5307 

-3.9669 

5.0609 

2.2915 

2.1000 

0.198 

2.7026 

1.5805 

-12.5921 

12.9758 

U.666U 

3.2000 

0.302 

-0.1828 

3.1257 

-38.3U07 

38.ii693 

3.0715 

li.2000 

0.396 

-2.7289 

1.5351 

-10ii.3707 

10U.U210 

-1.11*35 

5.3000 

0.1*99 

-2.6058 

-1.7358 

-313.609U 

313.6351 

-Iu3l59 

6.iiO00 

0.603 

0.3&7 

-3.1097 

-9i£.l558 

9ii2.1909 

-2.7099 

7.U000 

0.678 

2.8137 

-1.3729 

-2561.0508 

2561. 13U5 

1.5215 

8.5000 

0.801 

2.5000 

1.88U9 

-7693. 82U5 

769U.0709 

U.6313 

9.5000 

0.896 

-0.235ii 

3.1222 

-20913.981-3 

2091ii.6526 

3.5551 

10.5998 

1.000 

-2.8889 

1.2070 

-62816.530U 

62818.2167 

O.OOUli 
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Table  11.     Calculation  for  equation  (3-26),  0-  =   8.9133 


1 

2 

3 

h 

5 

6 

7 

k2£2 

x2/L 

A^inkgXg 

A2cosk2x2 

A-sinhk^- 

Ai  coshkpXp 

3*- Itf  5+6 

0.8913 

0.100 

-2.975it 

-0.628L 

-1.0138 

1.1*21*2 

-3.1931* 

1.7827 

0.200 

-3.7392 

0.2103 

2.8830 

3.0570 

-3.3S99 

2.7000 

0.303 

-1.63ii8 

o.90ia 

-7.1*02*1 

7.1735 

-O.6613 

3.6000 

O.iiOU 

1.6925 

0.8968 

-18.2800 

18.3128 

2.6221 

lu5ooo 

o.5o5 

3.7388 

0.2103 

-2*1*.9895 

1*5.012*1 

3.971*2 

5.^000 

0.595 

2.9559 

-0.63U7 

-110.6677 

110.7055 

2.3590 

6.2000 

0.696 

2.6128 

-0.9965 

-21*6.2996 

22*6.3755 

1.6922 

7.1000 

0.797 

-2.7881* 

-0.682*5 

-605.8013 

605.9839 

-3.2905 

8.0000 

0.898 

-3.732*2; 

0.1U55 

-12*90.0371 

11*90.2*791 

-3.1915 

8.9133 

1.000 

-I.8723 

0.8720 

-3713.9995 

3715.0001 

0.0002 
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Fifth  symmetrical  mode  —  0±      11.75U6,     02  =  9.8814 


Table  12.     Calculation  for  equation  (3-25),  0-,  =  U.751j6 


1 

2 

3 

I 

5 

6 

7 

h*i 

Xj/L 

AlSinklXl 

jUcosk-^Xn 

AoSinhk^x, 

ALCoshk-jX, 

3+U+5+6 

1.1755 

0.100 

1.0979 

-0.U581 

-X.7U3U 

2.1106 

I.012U 

2.1|000 

0.202 

O.8036 

0.8772 

-6.5026 

6.6105 

1.7887 

3.6000 

0.30k 

-0.526U 

1.0668 

-21.7521* 

21.781*9 

0.5729 

ii.7000 

0.397 

-1.1895 

0.01U8 

-65.39lit 

65.1*019 

-1.1639 

5.9000 

0.U97 

-0.U47 

-1.1031; 

-217.1227 

217.1259 

-l.'51*l*9 

7.1000 

0.598 

0.8671 

-0.8lli3 

-720.8775 

720.878U 

0.0537 

8.3000 

0.699 

1.0731 

0.5132 

-2393.3992 

2393.3991* 

1.5861 

9.5000 

0.801 

-0.0895 

I.I863 

-79U6.3655 

79k6.36$6 

1.0969 

10.6000 

0.892 

-1.0976 

0.!i586 

-23867.1i325 

23867. 1*325 

-0.6390 

11.75U6 

1.000 

-0.8631 

-0.S187 

-7571*0. 6186 

757U0.6186 

-0.01*1*1* 

5o 


Table  13.     Calculation  for  equation  (3-26),  02  =  9.88UU 

12  3  U  5  6  7 

k2x2            x2A  ^sink^  B2cosk2x2  B^inhk^  B^coshk^  3+U+5+6 

0.988U  0.100  3.5696  -C.5500  -1.1572  1.5296  3.3920 

1.9769  0.200  3.9270  0.3950  -3.5U05  3.679U  k.U609 

3.0000  0.303  0.6031  0.9903  -10.0169  10.0677  1.6i;39 

U.0000  O.i-Oii  -3.23ii9  0.6536  -27.2872  27.3082  -2.5603 

5.0000  0.505  -li.0988  -0.2837  -7lul958  7U.  2099  -6.9227 

5.9000  0.596  -1.5982  -0.9275  -182.U991  182.5201  -2.5QU7 

6.9000  0.697  2.U72U  -0.8l57  -I96.9873  196.1379  1.7073 

7.9000  0.798  li.2698  0.0U60  -73U8.5061  13I18.6U13  U.li^lO 

8.9000  0.899  2.litl5  0.865U  -3666.6201  3665.9868  3.3736, 

9.8810;  1.000  -1.8962  O.8962  -9809.9350  9810.9350  0.0000 
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First  anti-syxmetrical  mode  —  0-  -  1.6775,     02  _  l«kl06 

Table  Hi.     Calculation  for  equation  (3-27),     0X  =  1.6775'. 

12  3  k                    5  6  7 

kx              x/L  C  sink  x  C  cosk  x        C  sinhkx  C.  coshkx  3f  U*~5+6 

0.1678  0.100  -0.6922  3.3953  O.6989  -3.1*921  -0.0901 

0.3355  0.200  -1.36U6  3.25lii  I.U168  -3.6391  -0.3355 

0.5033  0.300  -1.9992  3.0162  2.1753  -3.8899  -0.6968 

0.6710  0.U00  -2.577U  2.6969  2.99l;3  -^.21*82  -1.1339 

0.8388  0.500  -3.0331  2.3013  3.8993  -U.7276  -1.6101 

1.0065  0.600  -3.5026  1.8U16  h.9132  -5.3U02  -2.0880 

1.1733  0.700  -3.8218  1.3330  6.0585  -6.0981*  -2.5287 

1.3lil0  0.800  -i;.0361i  O.7810i  7.3808  -7.0323  -2.9032 

1.5108  0.900  -U.1376  0.2066  8.9313  -8.18OO  -3.1791 

1.6775  1.000  -it.1215  0.3667          10.7055  -9.5368  -2.5861 


Uo 


Table  15.  Calculation  for  equation  (3-28),  0  =  1.U106 


1 

2 

3 

k 

5 

6 

7 

k2X2 

x2A 

D!sil*2X2 

D2cosk2x2 

D  sinhkgjCg 

D, coshk-x„ 
h          2  2 

3+U+5+6 

o.iiai 

0.100 

0.1651 

-0.9901 

-0.2330 

1.0100 

-0.0lt80 

0.2821 

0.200 

0.3269 

-0.96C5 

-O.I4.70U 

l.OUOl 

-0.0639 

0.U232 

0.300 

O.J482U 

-0.9118 

-0.7131 

1.0909 

-0.0558 

0.56l£ 

O.iiOO 

0.6280 

-0.3U50 

-0.9781* 

1.163U 

-0.0320 

0.7053 

0.500 

0.76l!i 

-0.761U 

-1.2591 

1.2592 

0.0001 

0.8i;6U 

0.600 

0.8796 

-0.6627 

-1.5650 

1.3801 

0.0320 

0.987U 

0.700 

0.9802 

-0.5509 

-1.9020 

1.528b 

0.0557 

1.1285 

0.800 

1.0615 

-O.I128O 

-2.2769 

1.7073 

0.0639 

1.2695 

0.900 

1.1216 

-0.2928 

-2.6971 

1.9200 

0.0517 

1.U106 

1.000 

1.1595 

-0.1597 

-3.1712 

2.1712 

0.0000 

la 


Second  anti-symmetrical  mode  —  <L-   li.7l8l,  £L  =  3.9680 


Table  16.  Calculation  for  equation  (3-27),  0,  =  Iw7l8l 


1 

2 

3 

h 

5 

6 

2 

klxl 

3^/L 

C  sink-x. 

CgCosl^x- 

C.sinhk.3c. 

Ci  coshk-x- 

3+1+5+6 

O.U179 

0.100 

1.51*53 

-3.0627 

-1.661.2 

3.8281* 

0.61*68 

0.9U37 

0.200 

2.7523 

-2.0177 

-3.7055 

5.0866 

2.1157 

l.lll56 

0.300 

3.3581 

-0.5316 

-6.5880 

7.U990 

3.7378 

1.8875 

0.1*00 

3.2303 

1.0707 

-10.981*8 

11.6122 

1*.  91*31* 

2.U000 

0.508 

2.6962 

2.5355 

-18.5807 

19.1071; 

.  5.7581* 

2.8000 

0.593 

1.1387 

3.2398 

-27.81*59 

28.3769 

li.9095 

3.3000 

0.698 

-0.5361 

3.3955 

-1*6. 0180 

h6. 6769 

3.5183 

3.8000 

0.801; 

-2.0800 

2.7199 

-75.9361 

76.8190 

1.59U8. 

li.3000 

0.910 

-3.11U3 

1.3782 

-125.2371 

126.7318 

-0.2U1U 

U.7187 

1.000 

-3.3992 

-0.0217 

-190.3783 

192.6100 

-1.1892 

k2 


Table  17.  Calculation  for  equation  (3-23),  0   =  3.9680 


1 

2 

3 

h 

5 

6 

7 

Xg/L 

D-,  sink  pXg 

DpCoskoXg 

DoSinhkpXg 

Di  coshlcpXp 

yhtfrs 

0.3968 

0.100 

-0.169U 

-0.9223 

-O.U203 

1.0798 

-O.ltfltf 

0.7936 

0.200 

-0.3125 

-0.7013 

-0.9136 

1.3318 

-0.5965 

1.1901; 

0.300 

-0.1*071 

-0.3713 

-1.5U79 

1.7962 

-0.5319 

1.5872 

O.l;00 

-O.U38U 

0.016U 

-2.1+332 

2.5173 

-0.3079 

1.98U0 

0.500 

-0.2iOl5 

0.1*015 

-3.6996 

3.70U6 

0.0050 

2.1i000 

0.605 

-0.2961 

0.7371* 

-5.6772 

$.^69 

0.3210 

2.8000 

0.706 

-0.1U69 

0.9ii22 

-8.5081 

8.2527 

0.5399 

3.2000 

0.807 

0.0256 

0.9983 

-12.7183 

12.2866 

0.5919 

3.6000 

0.908 

0.19U0 

O.8968 

-18.9912 

18.3128 

0.1a2U 

3.9679 

1.000 

0.322U 

0.6776 

-27.1^73 

26.10i62 

0.0011 

:,5 


Third  anti-symmetrical  mode  —  0.  =  7.3169,  02  =  6.1550 


Table  18.  Calculation  for  equation  (3-27),  fl^  =  7.3196 


1 

2 

3 

U 

5 

6 

7 

"ft 

XjA 

C  sink  x 

C2cosklXl 

C  sinhk  x 

C.  coshk  x 

3+U+5+6 

0.7320 

0.100 

r7.6l95 

8.U77ii 

9.109k 

-H;.5898 

-U.6225 

1.U639 

0.200 

-11.33U6 

1.2161; 

23.3202 

-25.9530 

-12.7513 

2.2000 

0.301 

-9.2166 

-6.707U 

50.8C92 

-52.0622 

-17.1770 

2.9000 

0.396 

-2.7268 

-11.0669 

103.2758 

-103.8827 

-lluUOOo 

3.7000 

0.506 

6.0395 

-9.6661 

230.i;007 

-230.6378 

-3.8637 

ij-.UOOO 

0.602 

10.61*79 

-3.5021* 

U6U.1837 

-U6U.23ii3 

7.29U7 

5.1000 

0.696 

10.5537 

U.3082 

93^.8573 

-93U.7U6U 

lii.9728 

5.9000 

0.806 

li.2623 

10.5711 

2080.625U 

-2080.251i6 

I5.20it2 

6.6000 

0.902 

-3.5510 

10.8298 

1*189.8876 

4il89.0950 

8.07lli 

7.3196 

1.000 

-9.8105 

5.80U7 

860li.lOi;7 

-8602.7609 

-2.3520 

Ui 


Table  19.     Calculation  for  equation  (3-28),     J0L  =   6.1550 


1 

2 

3 

k 

: 

6 

7 

k2*2 

XgA 

D^ink^ 

D2cosk2x2 

D^sinhkgXg 

Di  coshk  x„ 

3+ii+5+6 

0.6155 

oaoo 

-8.9983 

-0.8165 

-0.6579 

1.1955 

-9.2772 

1.2310 

0.200 

-Hi.  6927 

-0.3333 

-1.5730 

1.8583 

-1U.7I|07 

1.81.65 

0.300 

-2i;.9950 

0.2722 

-3.1032 

3.2U77 

-lii.5783 

2.5000 

0.U06 

-9.3271 

o.eon 

-6.0762 

6.1323 

-8.U699 

3.1000 

o.5oU 

-0.6U83 

0.9991 

-11.122a 

11.1215 

0.31*82 

3.7C00 

0.602 

8.2565 

0.3181 

-20.2982 

20.2360 

9.01*22* 

U.3C00 

0.698 

14.2782 

O.2;0O8 

-37.0015 

36.8567 

ll*.53**2 

U.9000 

0.796 

15.3111 

-0.1865 

-67.1299 

67.IU86 

IU.83U6 

5.5000 

0.89U 

10.992*6 

-0.7C87 

-122.8700 

122.3U08 

9.7639 

6.i55o 

1.000 

1.9916 

-0.9918 

-236.5ii52 

235.5Ui5 

0.0009 

is 


Fourth  anti-symmetrical  mode  ~  0^  =  8.291U,     02  =   6-9?22 


Table  20.     Calculation  for  equation  (3-27),     0.  =  8.2911; 


1 

2 

3 

1* 

5 

6 

7 

■ft 

^/L 

C  sink  x 

C  cosk  x 
2        11 

C  sinhk  x 
3          11 

C.coshk  x 
1*          11 

3+2*+5+6 

0.8291 

0.100 

-1.6925 

1.52*79 

2.1288 

-3.1288 

-1.11*28 

1.6583 

0.200 

-2.2868 

-0.20C5 

5. 8076 

-6.2101 

-2.9208 

2.5000 

0.301 

-1.3739 

-1.8379 

13.8882 

-ll*.0687 

-3.3923 

3.3000 

0.398 

0.3620 

-2.2655 

31.0762 

-31.11*33 

-1.9706 

it.  2000 

0.506 

2.0008 

-1.121*8 

76.5221 

-76.5132 

0.881*9 

5.0000 

0.602 

2.2012 

0.6509 

170.333U 

-170.252li 

2.9331 

5.8000 

0.698 

1.0665 

2.0315 

379.0979 

-378.8901 

3.3058 

6.7000 

0.796 

-0.9292 

2.0978 

932.1*371* 

-931.9121 

1.6936 

7.5000 

0.903 

-2.1532 

0.7952 

2075.1800 

-207U. 0061 

-0.181O 

8.2911; 

1.000 

-2.0788 

-0.9718 

1*518.8522 

-1*576.2591 

-0.1*575 

:,6 


Table  21.  Calculation  for  equation  (3-28),  0  =  6.9722 

1  2        3  h  5  6  7 

k2x2  x2^1  D-iSinkgXg  DpCoskpXp  D-sinhko^      Djcoshk^  3+1^5^6 

0.6972  0.100           1.7895  -0.7666  -0.7565  1.2530  1.51& 

1.39Hi  0.200           2.7U37  -0.1755  -1.8959  2.1U03  2.8126 

2.1000  0.301          2.I057  o.5oii8  -k.0295  h.ibk3  3.0253 

2.8000  0.it02            0.9336  0.91*22  -8.2075  8.2527  1.9210 

3.5000  0.502  -0.9776  0.9365  -16.57^0  16.5728  -0.01*23 

ii.2000  0.602  -2.1*291  0.1i903  -33.3990  33.3507  -1.9871, 

it.9000  0.702  -2.7381  -0.1865  -67.2688  67 .11*86  -3.0^8 

5.6000  0.802  -1.7591;  -0.7756  -135.1*683  135.2151  12.7882 

6.3000  0.903            0.01*65  -0.9999  -272.8023  272.2869  -1.1*685 

6.9722  1.000  - 1.7719  -0.7719  -53U.2861  533.2838  -C.C023 


hi 


Fifth  anti-syraietrical  node  —  01  -  10.9617,     02  =  9.2177 

Table  22.     Calculation  for  equation  (3-27),     0±  =  10.9617 

12,3  k  5  6  7 

k^           x^/L        C-jSink^  CgCosk^  C^inhk^  C^coshlc^  3+U+5+6 

1.0962  0.100            2.85ia  -l.i|66u  -li.2653  5-3378  2.1*602 

2.2000  0.201            2.$9h2  1.888U  -lli.30l5  Ui.Witf  U.6I46O. 

3.3000  0.301          -0.5060  3.1687  -U3.U391  1*3.5588  2.7801; 

l-.Uooo  o.Uoi        -3.0531;  0.9861  -130.6560  130.6996  -2.0237 

5.5000  0.501          -2.2637  -2.27lil  -392.561i9  392.5903  -!u5l2h 

6.6000  0.601            0.9995  -3.CU90  -1179.3890  1179.3890  -2.0082 

7.7C0O  0.707            3.1708  -0.1*922  -35U2.9625  35U3.07U2  2.7903 

8.8000  0.802           I.8768  2.6027  -105ii3.61i90  106U3.9S1U  li.8119 

9.9000  O.902          -1.U680  2.8533  -31975.2898  31976.2863  2.3818 

10.9617  1.000          -3.2068  1.0868  -92iii;9.5778  92l;52.1;590  0.7612 
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Table  23.  Calculation  for  equation  (3-28),  02  -  9.2177 


1 

2 

3 

k 

5 

6 

7 

k2X2 

VL 

J^sink^ 

D2cosk?x2 

D  sinhkgXg 

D,  coshlCpXp 

3+ii+5+6 

0.9218 

0.100 

0.0823 

-0.60l4i 

-1.0582 

1.U558 

-0.12UO 

1.8435 

0.200 

0.1001 

0.2693 

-3.0808 

3.238U 

0.5270 

2.8000 

0.302. 

0.03l;8    • 

0.9U22 

-3.1935 

8.2527 

1.0363 

3.7000 

0.U02 

-0.0551 

0.8U81 

-20.2153 

20.2360 

0.8137 

U.6000 

Oj;98 

-0.1033 

0.1122 

-^9.72*70 

U9.7U72 

0.0092 

5.5000 

0.596 

-0.0733 

-0.7087 

-122.363U 

122.3U80 

-0.802li 

6.5000 

0.70U 

0.022U 

-0.9766 

-322.6366 

322.5716 

-1.0192 

7.U000 

0.8d!t 

0.093k 

-0.it385 

-818.1555 

817.9925 

-0.5081 

8.3000 

0.902 

0.0937 

0.U31U 

-2012.3385 

2011.9363 

0,1232 

9.2177 

1.000 

o.ozlU 

0.9786 

-5037. 9iOU 

5036.9300 

0.003U 
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USE  RAYLEIOH'S  METHOD  TO  SOLVE  THE  NATURAL  FREQUENCY  OF  THE 
FIRST  NORMAL  MODE  OF  VIBRATION  OF  THE  PREVIOUS  FRAME 


Derivation  of  formula 


55 


L, 


El, 
(AS), 


#07? 


77777 


#7}77 


Fig.  17 

The  above  discussion  shows  that  the  first  normal  mode  of  vibration  is 
the  first  anti-symmetrical  mode  as  shown  in  Fig,  17.  Because  of  symmetry, 
the  shear  force  V  and  the  bending  moment  M  at  the  end  of  the  vertical  bars 
are  the  same. 
Determination  of  the  deflection  equation 

For  the  vertical  bar 

by  equation  (5) 


EIl-^-EIl^Il(xl)-V(Ll-V-M 


(U-l) 


integrate  equation  (U-l) 


EI, 
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A-  X „(xj  =    V(Lr     -  ix?)  -  Ms,  +  XL  (U-2) 


[id|xi(xi)=  7(llxi  -  *5)  ■  ^1  +  *1 

integrate  equation  (it— 2) 


2         3 


boundary  condition 

(1)  X^)       =  0  1^=0 

(2)  X^)       =0  n2  =  0 

^  -  0 

and  equation  (ii-2)  takes  the  form 

2 
E^C:^)  =  VC^  -  J.)  -  tt^  (U-ii) 

equation  (i:.-3)  takes  the  form 

T  2    3    ,,2 
L.X,   x,    1-cc, 

For  the  horizontal  bar 

by  eouation  (5)  » 

xo 
SI2X^(x2)  -  -  M  +  2M  -±  (h-6) 

integrate  equation  (li-6) 


■61 


x. 


EI2X£(x2)  =  -  I!x2+    M  —■  +  i^ 
integrate  equation  (l;-7) 


(ii-7) 


EI2Z2(X2)  =  -  —  +3lf  ^n3X2+nl 


thrB) 


n^,  n,  ,  are  integrating  constants  to  be  determined  by  boundary 

conditions 

Boundary  condition 


a)  x2u2J 


nu=  0 


::J 


=  0 


(2)  S2(x2) 


ML, 


=  0 


n_  = 


x2  =  L/2 


"3"  "5" 


with,  lu,  n»  .  thus  found,  equations  (U-7)  and  (h-&)  takes  the  form 


x        ML 
EE2X2(X2)=-I&2  +  Mlf +  ^ 


<1*~9) 


2  3 

I&2       Mxg       ML2 
EI2X2(x2)  =  _-2~+~-+-g~x2 


CU-io) 


canpatability  equation:     the  joint  at  b  is  rigid 


1M 


=  Z^(x2) 


"l"1! 


,x2  =0 


EEj- (~r  -  ml^  -^  (-^) 


simplify  and  rearrange 

1 


M  = 


(2  + 


VIX  =  PVL, 


1  21 
3  LXI2 


where 
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CHl) 


~  =  o  +  _  _ 


1  L2I1 


3  ¥2 


(!t-12> 


reflection  equation  for  vertical  bar 

Substitute  equation  (1;-12)  into  equation  (U-£)  and  rearrange 

2   ..3 


W-i[ft-«^->] 


Let  X_   be  the  value  of  X-.C2C.)  atx.  =  L 


vi: 


X1X1  "  E^ 


1  2-3P 


rearrange 


Oi-13) 


v_  fm,  6 

EI  "  *73  2^3P 

1 

substitute  equation  (U-lli)  into  equation  (it -13)  and  simplify 


CU-iM 
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W  -  ^[3(1  -  P)(^)2  -  £)3]  0.-15) 

equation  (U-15)  is  the  deflection  equation  for  vertical  bar 
differentiate  equation  (it— 1^ )  twice  and  simplify 


Deflection  equation  for  horizontal  bar 

Substitute  equation  (h-12)  into  equation  (li-10)  and  rearrange 

h<*t>-*m  (|)2  <sb>  [§>  - 1  <i|)2+  I  <J)3  £|>]  (U-17) 

differentiate  equation  (li-17)  tiri.ce  and  simplify 

L2       L2X1  ^ 

Determination  of  the  frequency  equation 

(1)  The  kinetic  energy  (K.S.)  and  the  potential  energy  (P.E.)  of  the 

vertical  bar 

\(\)  -  ^g  [3(1-P)(^)2  -  (%3]  (ii-15) 

r- 

X2 

xl<xl>  -  ^2  h  [9d-^2A  -  6(1-?)^  +  4]  (U-16) 
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Ii 

_ 


E.E.  =  |(A0)f2  /   x2^)  d^  (U-17) 


0 
substitute  equation  (ii-16)  to  equation  (it-17)  and  simplify 


K.E.  o  \   (Ap)Vx2__  _L^  (13  .  13  p   9  P2)  (u_l8) 

d  1   ill  (2_3P)2  15   5     5 


square  equation  (li-l6) 


^2   X2 

— - 


[x^) J  -  J^i  G^-)2  [L2(l-P)2  -  21^1-P)^  +  x2]      (fc-20) 

EI    L       ,2 
p.E.  , —±  J   L^Cxjl)  I  ^  (U-21) 

substitute  equation  (i;-20)  into  equation  (it-21 )  and  simplify 

2 

P.E.-  2k  J^.  ■■  36  0  (i  -  P  +  P2)  (U-22) 

2    1°   (2-3?)2  3 

(2)  The  kinetic  energy  (£.E.)  and  the  potential  energy  (P.E.)  of  the 
horizontal  bar 

1      2  "^2 
K.E.1=  ~  (Ad)2cd  T  I2(x2)   dx2  (li-23) 

substitute  equation  (i;-17)  into  equation  (ii-23)  and  simplify.  Note  here 
K.E.^  is  due  to  the  transverse  vibration  of  the  horizontal  bar 


a 


12  V  J"; 

K.B.2  »|  (Ao)^2^  (U-25) 

K»E«  d  K»S._  +  K»2«„ 

*'!  ) 

where  K.E^  is  due  to  the  longitudinal  vibration  of  the  horizontal  bar 

P.S.  =  |  EI2  .fL2[xj(x2)]2  d>:2  (h-26) 

substitute  equation  (i;-l8)  into  equa  ion  (!i-26)  and  simplify 

P.E.  =  \  EI  35L1   36p    1  (U-27) 

*   2   Lj     (2-3P)2  3 

Total  K.E.  =  equation  (U-18)  +  equation  (L-25) 


2®  zlLlLl 


-  !    2„2    -    f    2(Ao)l  ,33      13  _x  9  ,,2. 

Jl    2  ('3?~   T"  p      3"  P  > 

L(2-3?)       Ji?       *  * 


L                           iA  I2           2 
+  -1     (AP)0(1+      2_1  P 1__)1  (U-28)- 

IO.  2         T£l2  (2-3P)2  210 -! 

12 
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Total  P.E.  =  equation  (li-22)  +  equation  (ii-27) 


fl  X1L1   36    \2  ,n  A  1  V^  D2- 

2   l|  (2-3?)2 


[|-2P+(2  +  i^)p2] 


(U-29) 


Rayleigh's  ecraation    Total  X.2.  =  Total  P.S. 

Equating  equation  (h— 28)  and  equation  (U-29)  and  simplify 

H  |-2P+(2+i^)P2 

2  1  36  3  3  L^ 


(Ap)lLl  (2"3P)  2        ,33      13„      9  „2N      L2h  ,,       ^l  P2         . 

§■  %  -  *-  ? +  ^  p  )  +  l~i~  (1  +  ~n  — : — ~2> 

(2-3P)  12  I^Ig  210(2-3?) 

(U-30) 

Numerical  illustrative  example 

This  is  the  same  problem  as  that  on  page  16.  With 

L  L  (Ao) 

~  .  0.25,    ~  =  1,     and    _1  =  0.5 

x2  L2  (Ap)2 

substitute  into  equation  (ii-12)  to  find  P 

P=   0.i<8 
With  these  values,  substitute  into  equation  (h-30),  to  find  id  y 

2  EI1 

«>    =  8.8312 


<A?>1L1 


SI 

«>  =  2.9717   / i~p 

(Ap)^ 
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From  Table  1,  for  the  exact  solution,  we  have 

rsr 


CO  =  2.81UO 


(Ap)^ 


Hence,  the  Hayleigh's  method  is  xri.thin  an  error  of  5*6%.    For  most  practical 
cases,  this  error  is  insignificant  on  account  of  the  safety  factor  usually 
used  in  any  structural  design  and  Rayleigh's  method  is  a  good  approximation. 


6U 


APPENDIX 


F.     TSUN 
]  "AD    12.    PrilA,    DcLP.    I 

12    F  (2F20.8*     I  10) 

DC    21     1=1,    N 

fpHIA 
CA=CCSF(PHIA ) 
SA  =  SINF(PI 
"   ,=CCSF(PH  : 

I 

-XPF(- 

=  f: 

S  H  A  =  ( 

CHI   -  "■-  '-' 

.,  v> 
D1=(PH1A* ( ] .      - 

D?=  (Phil  A) )-( " 

X1=D1/D2 
\'l  =  {?,  .■■/(?. ' 


X?=(2.v.  *#C  .25  ) 
Y  2  =  (  C  B  - 1  .  ■:  )/r 
PU^CH    2 

fcr: 


21 


'  .  ?5  )  )  *  (  C!  *  C    •*.  -1  i  0  )  /  ( '  ~     '  *CHA  ) 

■  )  /  (8.C*Si 
[  1  .(> 

lIBf     Yl  »     Y2,     XI.     X2 

a  t    (  3  F  2  n .  .  a  ; 

PHI A=P  I  . 

STOP 

END 
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ABSTRACT 

By  applying  the  Bernoulli-Euler  equation,  the  equation  of  motion  for  a 
vibration  elastic  beam  is  derived  in  the  beginning  of  this  report.  Since  a 
rectangular  rigid  frame  can  be  considered  as  three  elastic  beams  joined 
rigidly  at  right  angles,  the  frequency  equation  can  be  obtained  from  direct 
expansion  of  the  frequency  determinant  that  follows  as  a  consequence  of 
satisfying  boundary  conditions  of  the  equations  of  motion  for  the  vibrating 
elastic  beams.  In  order  to  show  how  to  use  these  frequency  equations,  a 
numerical  example  is  presented  and  the  natural  frequencies  and  mode  shapes  of 
the  first  ten  modes  of  vibration  of  a  simple  rectangular  fixed  end  frame  are 
obtained  by  the  aid  of  a  1620  digital  computer. 

Rayleigh's  method  is  also  used  to  find  the  natural  frequency  of  the 
first  vibration  mode  for  the  same  numerical  example  and  the  error  compared  to 
the  exact  solution  is  found  to  be  5*6£. 


